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Abstract. Cold, velocity-controlled molecular beams consisting of a single quantum 
state promise to be a powerful tool for exploring molecular scattering interactions. 
In recent years, Stark deceleration has emerged as one of the main methods for 
producing velocity-controlled molecular beams. However, Stark deceleration is shown 
not to be effective at producing a molecular beam consisting of a single quantum 
state in many circumstances. Therefore, quantum state purity must be carefully 
considered when using Stark decelerated beams, particularly in collision experiments 
where contributions from all quantum states must be addressed. 



Rotational- state purity of Stark- decelerated molecular beams 



2 



1. Introduction 

Molecular beams have been the workhorse for studying gas-phase interactions since 
the 1950's [D l2]- A recent goal of these studies has been to measure and possibly 
even control interactions on a state-by-state level, allowing exploration of the role of 
quantum mechanics in molecular processes. A prerequisite for these studies is the 
ability to prepare the reactants with well defined internal (electronic, rovibrational) 
and external (translational) degrees of freedom. Examples of experiments that exploit 
restricted degrees of freedom include controlled chemistry [3llll[5], stereochemistry [6l[7], 
fundamental studies of bimolecular reaction dynamics |8l |9], and scattering resonances 
[T0| [TT] . In particular, studies of scattering resonances benefit enormously from the 
use of beams with narrow velocity distributions [121 [13] where the collision energy can 
be varied over the resonance [H]. The study of cold molecules [I5] is predicted to be 
particularly insightful for exploring intermolecular interactions and chemistry [T6| [T7]. 
For example, increased control of the participating quantum states can elucidate the role 
of barriers in molecular reactions where ground-state interactions may be energetically 
forbidden but are possible among excited states. 

Several methods have been developed to produce cold molecular samples with 
limited internal and external degrees of freedom. These methods include supersonic 
expansion [18], buffer-gas cooling [IHIEO], velocity selection [21], kinematic cooling [22], 
photoassociation [231 121], Feschbach resonances [251 126], and beam deceleration using 
optical [27], magnetic [28] or electric fields [29]. One such beam deceleration technology. 
Stark deceleration, uses inhomogeneous time- varying electric fields to decelerate a beam 
of polar molecules produced via supersonic expansion. This method has been used 
to decelerate a variety of molecular species, including CO, NH3, ND3, OH, SO2, NH, 
and CaF [291 EH EH 1321 [331 [31] and has been employed in a number of cold collision 
experiments [H ESI ESI EH EHl [33 [IQ] . 

There are numerous experiments that require samples to be both translationally 
cold and in a well-defined single internal quantum state. Having proven to be highly 
effective at producing cold molecules. Stark deceleration is a logical candidate for the 
production of these cold molecular beams consisting of a single state. Therefore, 
the physical chemistry community would benefit from a systematic study specifically 
addressing the state selectivity of Stark deceleration. 

Explorations of Stark deceleration selectivity have appeared previously in the 
literature. One study focused on the deceleration of different isotopologues [30] while 
another involved different projections of the dipole moment of the OH molecule onto 
the electric field axis ^Ij. Both of these studies explored selectivity among species 
with different Stark shift to mass ratios, which is the parameter that determines the 
dynamics of the deceleration process. The present investigation expands on this previous 
work with a rigorous and detailed experimental study of the state selectivity of Stark 
deceleration. The principal method used for this investigation is the deceleration and 
detection of various rotational states of a single isotopologue of ammonia. Specifically, 
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molecules in two rotational states, which have different Stark shifts, are decelerated and 
state-selectively detected using decelerator timing sequences corresponding to a large 
range of Stark shift to mass ratios. These measurements are compared to the results 
of a series of three-dimensional Monte Carlo simulations to confirm the validity of the 
computational model. This model is then used to extract the final mean velocity for 
different rotational states decelerated under various conditions. We find that for a given 
timing sequence, the final mean velocity varies with rotational state. This velocity is an 
important parameter for many collision experiments where collision energies must be 
precisely characterized. 

2. Stark Deceleration 

The Stark deceleration process begins with a pulsed supersonic expansion, where 
the molecules of interest are seeded as a low concentration in a heavier buffer gas. 
The resulting molecular beam is internally cold with molecules predominantly in the 
rovibrational ground state. However, if the molecule has multiple low lying, energetically 
accessible excited rotational states, as is the case with ammonia, a significant portion 
of the population can end up in these excited states. For molecules such as OH, where 
the first excited rotational state is ~ 80 cm~^ above the ground state, only the ground 
state will be substantially populated and rotational excitations can usually be neglected. 
Though internally cold, such beams have a high forward velocity in the laboratory frame 
of many hundreds of meters per second, which can be reduced by Stark deceleration. 

The Stark deceleration technique [29] relies on the interaction of an electric field 
with a molecule's electric dipole moment to reduce the high forward velocity a portion 
of the molecular beam. A series of high-voltage electrode pairs create periodic electric 
field maxima along the beam axis. Any molecule whose energy increases with increasing 
electric field (known as a weak-field-seeking state) will experience an acceleration 
opposing its motion. As these molecules propagate into the increasing electric field, 
longitudinal kinetic energy is converted into potential energy. If the molecules were 
allowed to continue past the local electric field maximum they would regain their lost 
kinetic energy. However, before they reach the field maximum, the electric field is 
switched off rapidly, reducing the forward velocity of the beam by a small amount. This 
process is repeated using successive electrode pairs to produce a molecular beam with 
the desired final velocity. Transverse guiding is also achieved because the molecules 
are attracted to the electric field minimum along the center of the decelerator. To 
guide equally in both transverse dimensions, successive electrode pairs are oriented 
orthogonally to one another. 

The final velocity of the decelerated molecules is determined by the amount of 
energy removed per stage and the total number of stages. The amount of energy removed 
per stage can be described by the position of the so-called "synchronous molecule" , which 
is at the same physical location relative to each electrode pair at the time the electric 
fields are switched off. This position is referred to as the phase angle, 0o, where 0o = 0° 
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Figure 1. Rotational energy levels of the ground electronic state in ^^NDs. The states 
are labeled by \J,K) and organized by K value. The inversion splitting is not visible 
on this scale. 



corresponds to no energy being removed and (f)Q = 90° corresponds to the maximum 
amount of energy being removed. Increasing the phase angle, and hence the amount 
of energy removed per stage, results in both lower final velocities and lower overall 
molecule numbers [12]. This is because the phase-space acceptance, the region of phase 
space that will be decelerated and transported to the end of the decelerator, decreases 
with increasing phase angles. 

Molecules close to the synchronous molecule in phase space experience an average 
force toward its location. This force results in an effective time-averaged moving 
potential well [SO]. The dynamics of molecular trajectories in this moving potential 
well have been thoroughly investigated in the case where the molecules have the same 
mass and Stark shift as that of the synchronous molecule |13]. In the present analysis, 
molecules with the same mass but different Stark shifts are experimentally investigated. 
This required a molecular beam with many initially populated weak-field-seeking states. 
The numerous low-lying rotational states of the ammonia molecule (Figure [1]) are ideal 
for this purpose because of their appreciable Boltzmann factors at our beam temperature 
(6-7K). 

3. The Ammonia molecule 

In its rovibrational ground state, ammonia has a trigonal pyramidal structure with a 
ladder of rotational states denoted by \J,K), where J is the total angular momentum and 
K is the projection of J onto the symmetry axis of the molecule. The trigonal pyramidal 
structure has two inversion eigenstates, which are symmetric and antisymmetric 
combinations of the two polarized "umbrella" states. The umbrella states are defined 
by the nitrogen being either "above" or "below" the plane defined by the hydrogen 
atoms. Because there is a barrier for the nitrogen to tunnel through this plane, the two 
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Table 1. The energies, inversion splittings, and relative dipole moments for various 
rotational levels. Entries in the fourth column include all possible non-zero values of 
M. 



\J,K) E,ot (cm-i) AU„, (cm-i)[ 



1 1 

2 2 

2 1 

3 3 
3 1 
3 2 



8.26 
22.74 
28.82 
43.51 
53.60 
59.69 



0.0530 
0.0525 
0.0520 
0.0518 
0.0511 
0.0506 



Me// 



1/2 
2/3, 1/3 
1/3, 1/6 
3/4, 1/2, 1/4 
1/4, 1/6, 1/12 
1/2, 1/3, 1/6 



Ionization 
laser 



MicroChannel Plate 



Pulsed 
valve 



Decelerator 
electrodes 



Skimmer 



Mass 
spectrometer 



Figure 2. Experimental set-up consisting of a PZT-driven pulsed valve, molecular 
beam skimmer, decelerator, time-of-flight mass spectrometer, and microchannel-plate 
ion detector. The decelerator consists of 150 electrode pairs (not all are shown). 



eigenstates are separated in energy by an inversion splitting, At/j^t,, which varies only 
slightly among low rotational states (Tabled]). When an electric field is applied, the two 
opposite-parity inversion states mix and are repelled from one another in energy. This 
Stark energy is given by 



stark 




E 



MK 

j(7TT 



(1) 



where /i is the electric dipole moment {fiND-i = 1-47 Debye), E is the electric 
field, M is the projection of J onto the electric-field axis, and the upper (lower) signs 
correspond to weak- (strong-) field-seeking molecular states. Fully deuterated ammonia 
(ND3) was chosen over NH3 because of its favorable Stark shift due to its much smaller 
inversion splitting (~ 0.05 cm~^ compared to ~ 0.79 cm~^ for NH3). 



4. Experiment 

The experimental apparatus consists of a pulsed valve. Stark decelerator, and molecular 
detector comprised of a pulsed dye laser, time-of-fiight mass spectrometer (TOFMS) 
and microchannel plate (MCP) (Figure E]). A molecular beam is created by supersonic 
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expansion of ^^NDs seeded in krypton via a PZT-actuated pulsed valve. The valve 
nozzle is constructed from a 2 cm long piece of hypodermic needle tubing and is 
resistively heated to approximately 185° C to increase the population of molecules in 
excited rotational states. The molecular beam, with an initial mean velocity of 460 m/s, 
passes through a skimmer into a differentially pumped vacuum chamber containing the 
149-stage Stark decelerator operating at ±12kV. After deceleration, the molecules are 
detected via a |J,K) state sensitive 2 + 1 resonance enhanced multi-photon ionization 
(REMPI) scheme near 317 nm. The ions are then accelerated by the TOFMS onto 
the MCP. The signal at the arrival time for an ion with a mass-to-charge ratio of 20 is 
measured and averaged over many experimental runs. 

Operation of the Stark decelerator begins by setting the high-voltage switching 
timings appropriate for the desired final velocity. This timing sequence is calculated 
from trajectory simulations using electric fields generated from Comsol, a finite element 
analysis package. To facihtate the discussion of how these timing sequences are related 
to different rotational states, we introduce the effective dipole moment, 

Mf = IH J^y (2) 

For the | J,K,M) = |1,1,1), |2,2,2), and |3,3,3) states, ^ = i |, and |, respectively. 
Although fieff does not uniquely determine the rotational state of the molecule, it does 
uniquely determine the timing sequence. Therefore, to slow synchronous molecules 
with different fieff values to the same final velocity requires operating the decelerator 
using different timings. For example, with our particular apparatus, decelerating a 
|1,1,1) synchronous molecule (^^^ = i, see Table [I]) from 460 m/s to 60 m/s requires 
calculating decelerator timings for a phase angle of 0o ~ 84°. Removing the same 
amount of energy per stage for a synchronous molecule with = 5 requires a much 
smaller phase angle of 0o ~ 7°. Calculated switch timings for decelerating the |2,2,2) 
and |3,3,3) states to 60 m/s, relative to the |1,1,1) state, appear in Figure |3l Even 
though the effective dipole moments vary from those of the 1 1,1,1) state by as much 
as 50 %, the resulting switch timings differ by at most 2.5 %. The inset in Figure [3] 
contains the relative timings for a hypothetical ND3 molecular state with = 5. 

Note that while the Stark shift depends on M, our molecular detection scheme is 
M-independent. Hence, when discussing timings for the decelerator, the full |J,K,M) 
of the synchronous molecule used to calculate the timings is referenced. Discussion 
regarding the actual detection of molecules references only |J,K), as all possible values 
of M are included. 

When evaluating the state selectivity of the Stark decelerator we consider only the 
relative populations in weak-field-seeking states. States that are strong-field-seeking 
as well as those unperturbed by the electric field will effectively be filtered out. The 
former are attracted to the highest electric fields, located at the surface of the electrodes, 
and therefore do not make it to the end of the decelerator. The latter are unguided 
transversely and not decelerated longitudinally and hence make up a portion of the 
undecelerated background molecular beam. 
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Figure 3. Percent differences for the ^^^^^ ~ | (|2,2,2), squares) and the | (|3,3,3), 
circles) timing sequences from the ^^^^^ ~ i (|1,1,1)) for decelerating to 60 m/s. Inset: 
timing difference of the — 5. 



To test how the final molecule density depends on the chosen timing sequence, 
we measure relative populations of decelerated rotational states using timing sequences 
for synchronous molecules with = |,|,|,2,3,4, and 5 decelerated to 60 m/s. A 
final velocity of 60 m/s is the lower limit for the 1 1,1,1) synchronous molecule given our 
beam's initial forward velocity and our number of decelerator stages. This constitutes 
a removal of more than 98% of the initial longitudinal kinetic energy of the molecular 
beam. This velocity was chosen for these measurements because a low final velocity 
results in decelerated molecules that are spatially displaced from the undecelerated 
background molecular beam. This displacement effectively amplifies any state selectivity 
of the decelerator as the undecelerated background molecular beam has a state purity 
determined solely by the dynamics of the supersonic expansion. 

5. Results and Discussion 

The distribution of rotational states initially present in the molecular beam was 
determined using the setup shown in Figure [21 but with no voltage applied to 
the decelerator (free fiight). Though not a thermal distribution, the approximate 
temperature was found to be 6—7 K by fitting the measured spectrum (Figure 111(A)) to a 
theoretical spectrum weighted by Maxwell-Boltzmann factors. The theoretical spectrum 
was derived from the Hamiltonian for ND3 using constants together with individual line 
strength calculations from |l5l |l6]. As shown in Figure HI (A), the initial molecular 
beam consists of numerous excited rotational states. The experimental spectrum has 
power broadened linewidths of ~ 0.25 cm~^, which is narrow enough that, for our chosen 
detection wavelengths, adjacent lines do not overlap, thus allowing for true single-state 
detection. Figure 111(B) shows the densities of molecules in states |1,1), |2,2), |2,1) and 
1 3,3) as a function of detection time under free fiight conditions. The traces have been 
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Figure 4. Characterizations of the initial molecular beam. (A) Measured REMPI 
trace of the molecular beam rotational state distribution where numerous rotational 
states are populated. The relative state populations indicate an approximate 
temperature of 6-7 K. (B) Time-of-flight signal of various rotational states at the end 
of the decelerator where no voltage is applied, as detected on the starred wavelengths 
in (A). Note that the |1,1) state (solid squares) has been divided by 10 for ease of 
comparison. 
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Figure 5. (A) Detecting |1,1) (squares) and |2,2) (circles) molecules with timings 
calculated for decelerating a ^^^j^ = ^ synchronous molecule to 60 m/s. (B) Same as 
A but using a ^^^^ = | timing sequence. (C) Same as A but using an = 5 timing 
sequence. Note the comparable ratios and the shift in arrival times. 



normalized to the calculated transition line strengths and are therefore on a common 
vertical scale to facilitate direct comparison. The |1,1) trace has also been divided by 
10 for display purposes. Relative rotational-state populations initially present in the 
molecular beam are quantified using ratios of peak molecule densities with the initial 
free-flight ratio of |1,1):|2,2) being ~ 22:1. 

The ratio of |1,1):|2,2) is also measured after the deceleration process. If Stark 
deceleration were capable of generating a molecular beam consisting of a single rotational 
state, this ratio would exhibit a strong dependence on the timing sequence used. For 
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instance, the ratio would be maximized when using timings calculated with = | 
(i.e. for a |1,1,1) synchronous molecule). Conversely, this ratio should approach zero 
when using timings determined with ^^^ll = | (i.e. for a 1 2,2,2) synchronous molecule). 
Observation of no dependence on the timing sequence would indicate no purification 
of these two states by the Stark deceleration process. Measured time-of-fiight molecule 
densities of the |1,1) and |2,2) states are shown in Figure |5] for timings calculated using 

= |, |, and 5 synchronous molecules slowed to 60 m/s. The ratio of the density 
of 1 1,1) to |2,2) molecules does not change appreciably, even when the timing sequence 
for a dipole moment an order of magnitude larger is used. Thus, deceleration has little 
effect on the relative populations of these weak-field-seeking states as no correlation 
between state purity and timing sequence is observed. 

Comparing the ratio of |1,1):|2,2) between free flight and deceleration using = | 
timing we see an increase of a factor of 2 (from 22:1 to 41:1). The same increase in the 
|1,1):|2,2) ratio is seen when using = | timing (from 22:1 to 40:1), which indicates 
the effect is not due to the state selectivity of the deceleration process. This increase 
is primarily due to all M sublevels being present in the free-flight beam, but not in the 
decelerated beam. Specifically, 1/5 of the |2,2) molecules (|2,2,2) only) are decelerated 
compared to 1/3 of the |1,1) molecules (|1,1,1) only). The weak-field-seeking |2,2,1) 
molecules are absent in the decelerated beam because their Stark shift is too small to be 
decelerated down to these velocities. In addition to the M sub-level contribution, the 
|2,2) molecules have a 50 % wider initial velocity distribution than the |1,1) molecules 
after the supersonic expansion, decreasing the efficiency of the deceleration process and 
increasing the ratio of |1,1):|2,2). 

To further explore the effect of various timing sequences on state selectivity, a single 
state was detected using timing sequences for = i, |, |, 2, 3, 4, and 5 (Figure |6]). 
Increasing molecular densities are observed when timing sequences for larger are 
used, regardless of the state being detected. This behavior can be compared to results 
from a series of three-dimensional Monte Carlo simulations. The simulated density as a 
function of time-of-fiight of |1,1) and |2,2) molecules is shown in Figure [6] (B) and (D), 
respectively. The simulations correctly predict both the shift in arrival times as well as 
the increase in molecule density. 

The simulations were also used to extract information about the final beam that is 
difficult to determine experimentally. Of particular interest are the final molecular 
velocity distributions at the exit of the decelerator under different timing schemes 
(Figure [7]). As shown, decelerator timing sequences based on increasing result 
in increasing final velocities. This shift in velocity, together with the timing difference 
of the last stage (Figure [3]), combine to produce the observed shift in the arrival times 
(Figure [6]). Figure [7] (B) shows the final mean longitudinal velocity for both |1,1) 
molecules and |2,2) molecules. The resulting final velocities for these two states are 
nearly identical under all timing conditions. In a single realization of the experiment, 
only one timing sequence can be used, resulting in a molecular beam consisting of both 
|1,1) and |2,2) molecules with nearly identical final velocities and spatial distributions. 
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Figure 6. Comparison of experiment and simulation. Time-of-fliglrt traces are labeled 
by the corresponding ^^^^ value used to calculate the decelerator timing sequence. (A) 
Detecting |1,1) using various timing sequences calculated using synchronous molecules 
with different //e// values. (B) Results from 3D Monte Carlo simulations of the 
expected |1,1) densities. (C) Detecting |2,2) under the same conditions as in A. (D) 
Results from simulations for |2,2). 
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Figure 7. Resulting final velocity distribution for molecules in the |1,1} state using 
various ^'"^^ -value timing sequences as produced by Monte Carlo simulations. (A) 
Curves are labeled by the corresponding value used to calculate the decelerator 
timing sequence. (B) Mean longitudinal velocities extracted from simulations. 



Thus both weak-field-seeking states will be present in the resulting decelerated molecular 
beam. 

Studying gas-phase molecular interactions where the reactants are in a single 
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internal quantum state and have well-controlled external degrees of freedom is a major 
goal in physical chemistry. Stark deceleration is a possible candidate for producing 
these single-state molecular beams as it is a proven method for controlling external 
degrees of freedom. In addition, Stark deceleration is capable of filtering out strong-field- 
seeking states, states that are unperturbed by the electric field, and weak-field-seeking 
states with Stark shifts too weak to be slowed to the desired final velocity, making it 
an increasingly viable candidate. Our results show that in many circumstances Stark 
deceleration is not capable of filtering out all but one of the remaining weak- field-seeking 
states. In fact, every weak-field-seeking state present in the supersonic expansion that 
can be decelerated to the desired final velocity, will be decelerated regardless of whether 
its Stark shift is stronger or weaker than that of the synchronous molecule. This is an 
important consideration when interpreting results of collision experiments using Stark 
deceleration and analogous techniques. 
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